Draft version July 1 , 2008 

Preprint typeset using I4TgX style emulateapj v. 10/09/06 



oo 
O 
O 
(N 



REVEALING THE HIGH-REDSHIFT STAR FORMATION RATE WITH GAMMA-RAY BURSTS 
Hasan Yuksel 1 - 2 , Matthew D. Kistler 1 - 2 , John F. Beacom 1 - 2 ' 3 , and Andrew M. Hopkins 4 

Draft version July 1, 2008 

ABSTRACT 

While the high-z frontier of star formation rate (SFR) studies has advanced rapidly, direct measurements 
beyond z ~ 4 remain difficult, as shown by significant disagreements among different results. Gamma-ray 
bursts, owing to their brightness and association with massive stars, offer hope of clarifying this situation, 
provided that the GRB rate can be properly related to the SFR. The Swift GRB data reveal an increasing 
evolution in the GRB rate relative to the SFR at intermediate z; taking this into account, we use the highest- 
z GRB data to make a new determination of the SFR at z = 4 — 7. Our results exceed the lowest direct 
SFR measurements, and imply that no steep drop exists in the SFR up to at least z ~ 6.5. We discuss the 
implications of our result for cosmic reionization, the efficiency of the universe in producing stellar-mass black 
holes, and "GRB feedback" in star-forming hosts. 
Subject headings: gamma-rays: bursts — galaxies: evolution — stars: formation 
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1. INTRODUCTION 

The history of star formation in the universe is of intense 
interest to many in astrophysics, and it is natural to pursue 
pushing the boundary of observations to as early of times 
as possible. Our understanding of this history is increas- 
ing, with a consistent picture now emerging up to redshift 
z ~ 4, as summarized in Fig. Q] The cosmic star for- 
mation rate (SFR) measurements from the compilation of 
iHopkins & Beacoml (2006) are shown, along with new high-z 
measurements based on obser vations of color-selected Lyman 
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Break Galaxies (LBG) dBouwens et al.ll2008llMannucci et al. 



20071: IVerma et al]|2007l) and Lva Emitters (LAE) (Ot aetal. 



2008). Much current interest is on this high-z frontier, where 



the primeval stars that may be responsible for reionizing the 
universe reside. Due to the difficulties of making and in- 
terpreting these measurements, different results disagree by 
more than their quoted uncertainties. 

Instead of inferring the formation rate of massive stars 
from their observed populations, one may directly mea- 
sure the SFR from their death rate, since their lives are 
short. While it is not yet possible to detect ordinary core- 
collapse supernovae at high z, long-duration gamma-ray 
bursts, which have been shown to be ass ociated with a spe- 
cial class of core -collapse supernovae dStanek et al.l 120031: 
iHjorth et all 120031) . have been detected to z = 6.3. The 
brightness of GRBs across a broad range of wavelengths 
makes them promising probes of the star formation his- 
tory (SFH) (see, e.g., the early works of rlo tani 1997; 
Wiiers et al. 1998; Lamb & Reichart 2000; Blain & Natarajan 
2000; Porciani & Madau 2001; Bromm & Loeb 2002). In the 
last few years, Swift 5 ( Gehrels et al.N2004l) has spearheaded 
the detection of GRBs over an unprecedentedly-wide redshift 
range, including many bursts at z > 4. Surprisingly, ex- 
amination of the Swift data reveals that GRB observations 
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FIG. 1 . — The cosmic star formation history. The compiled SFR data (light 
circles) and fit (dotted line) of Hopkins & Beacom 1 2006) are shown, as well 
as newer high-z data (the LAE points only sample Lya Emitters). The results 
of this work, as inferred using bright Swift gamma-ray bursts, are shown with 
dark diamonds. The solid line is our new high-z fit given by Eq.|5] 



are not tracing the SFH d irectly, instead implying some kind 
of additional evolution (Daigneetal. 2006; Le & Dermer 
2007; Yuksel & Kistler 2007; Sal vaterra & Chincarini 2007; 
Guetta&Piran 2007r iKistler et al.l 120081: ISalvaterra et all 
20081) . 



GRBs can still reveal the overall amount of star formation, 
provided that we know how the GRB rate couples to the SFR. 
In this study, we use the portion of the SFH that is sufficiently 
well-determined to probe the range beyond z ~ 4. We do this 
by relating the many bursts observed in z ~ 1 — 4 to the corre- 
sponding SFR measurements, and by taking into account the 
possibility of additional evolution of the GRB rate relative to 
the SFR. This calibration eliminates the need for prior knowl- 
edge of the absolute conversion factor between the SFR and 
the GRB rate and allows us to properly relate the GRB counts 
at z ~ 4 — 7 to the SFR in that range. Additionally, we make 
use of the estimated GRB luminosities to exclude faint low-z 
GRBs that would not be visible in our high-z sample, i.e., to 
compare "apples to apples". 
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Our results show that the SFR must be relatively high in 
the range z = 4 — 7 when compared to SFR measurements 
made using more conventional techniques. While the GRB 
statistics at high z are relatively low, they are high enough, 
permitting an approach complementary to other SFR deter- 
minations, which themselves must contend with presently un- 
known extinction corrections, cosmic variance, and selection 
effects, most importantly that flux-limited surveys necessarily 
probe the brightest galaxies, which may only contain a small 
fraction of the star formation activity at early epochs. We 
discuss the implications these results in the context of reion- 
ization and broader applications. 

2. THE GRB TECHNIQUE 

The relationship between the comoving GRB and star for- 
mation rate densities can be parametrized as nQRB(z) = 
£(z) x p*(z), where £(z) reflects the fraction of stars that 
produce long-duration GRBs and any additional evolutionary 
effects. Import antly, the SFH is well known for z < 4, and 
we use the Hopkins & B eacoml d2006l) fit in this range. The 
expected (all-sky) redshift distribution of GRBs can be cast as 



dN 

dz 



F{z 



£{z) p*(z) dV/dz 

(/beam) 1 + Z 



(1) 



where < F(z) < 1 summarizes the ability both to de- 
tect the initial burst of gamma rays and to obtain a redshift 
from the optical afterglow, and includes observational fac- 
tors such as instrumental sensitivities. Beaming rendering 
some fraction of GR Bs unobservable from Earth is accounted 
for through ( f h „ am ) (fBloom et alj|2003t iFirmani eld] 12004 
iKocevski & But ler 2008), the l/(l+z) is due to cosmological 
time dilation of the observed rate, and dV/dz is the comoving 
volume per unit redshift 6 . 

Both F(z) and £{z) are discussed in detail in lKistler etal.1 
(2008), in which it has been shown that the former can be 
set to a constant (Fq) by focusing on a bright subset of 
the GRB sample, allowing the latter to be parametrized as 
£(z) — £q(1 + Z ) Q ' w i m a — 1-5 ( an enhanced evolu- 
tion of GRBs) preferred over a = 0. Here £o is a (un- 
known) constant that includes the absolute conversion from 
the SFR to the GRB rate in a given GRB luminosity range. 
The evolutionary trend described by a may arise from sev- 
eral mechanisms dKistler et al.l 120081) . su ch as a GRB pref- 
erence for low-metallicit y e nvironm ents (Staneket al. 2006; 
lLanger & Normanl 120061: O l2007t ICen & Fand 12007b . and 
must be accounted for to characterize nGRB- 

Fig. |2] shows the luminosity-redshift distribution of the 
63 lon g-duration (Tgo > 2 s) GRBs from the catalog of 
iButler et al] J2007I) . We compute source-frame GRB lumi- 
nosities, L[ SO = _Ei so /[Tgo/(l + z)], based on the values 
of .Biso, the isotropic equivalent (uncorrected for beaming) 
1 — 10 4 keV energy release in the rest frame, and Tgo, the 
interval containing 90% of the prompt GRB emission. Using 
another quantity, such as the peak isotropic equivalent flux 
(N. Butler, private communication) yields quantitatively sim- 
ilar results. 

As only very bright bursts can be seen from all redshifts, we 
define a cut, based on the luminosities of the high-z events, of 
L[ so > 3 x 10 51 erg s _1 , for the subsample that we use to 



6 dV/dz = 4-n-(c/H )dl(z)/y/(l + z) 3 n m + C A , where d c is the 
comoving distance, Q m = 0.3, Qa = 0.7, and Hq = 70 km/s/Mpc. 



52 



io 53 

j- 10 



J 50 

10 



10 



49 



T 



1 J, 1 I 
* lit 



( 





■ Absorption 

Emission 
» Photometric — 

ii I i i i i I i i i 







1 



4 
z 



7 



FIG. 2. — The luminosity-redshift distribution of 63 Swift GRBs (two z < 
0.1 GRBs are not visible), as we determine from the data of Butler et al. 
(2007), marked according to the primary method of redshift determination. 
The ratios of the number counts of GRBs in the three shaded boxes are used to 
estimate the SFR density at high z (4-5, 5-7) by normalizing to the z = 1 — 4 
SFR, as discussed in the text. The boxes contain (21,4,4) bursts, respectively. 

estimate the high-z SFR 7 . (A naive normalization using all 
GRBs in z = — 4 without regard to luminosity would yield 
values ~ 2 times smaller.) This circumvents the need for an 
exact definition of the z-dependent threshold for detecting a 
GRB with a redshift measurement. The complicated detection 
threshold of Swift (e.g., Band 2006) and the human factor in- 
volved with optical observations makes a detailed treatment 
based on constructing the GRB luminosity function at low 
z and marginalizing over this distribution to obtain the SFR 
at high-z practically challenging, if not impossible. The net 
effect would probably only be to make our eventual results 
somewhat larger. 

The shaded boxes in Fig.|2]show three groups of GRBs de- 
fined by this I/jso cut in z = 1 — 4, 4 — 5, and 5 — 7 (the high-z 
boxes chosen to have equal counts). The GRBs in z = 1 — 4 
act as a "control group" to base the GRB to SFR conversion, 
since this range has both good SFR measurements and GRB 
counting statistics. We calculate the "expected" number of 
GRBs in this range as 



.A/n = A^ f dzF(z)£(z 



A / dzp\(z){l + z) 



p*(z) dV/dz 

(/beam) 1 + 2 

a dV/dz 
1 + z ' 



(2) 



in which A = At Ail £q i 7 b/47r(/b ca m), an unknown normal- 
ization, depends on the total live-time, At, and the angular sky 
coverage, Ail. Using the average SFR density, (p*) z _ z , a 
similar relation can be written for z = 4 — 5 and 5 — 7 as 

m::? Z2 = (p*) zi _ Z2 a f Z2 dz (i + zr ^ . o) 

J Z\ J- + Z 

Our interest is in {p*) z _ z , which we find by dividing out A, 
using Eq.[3] Taking the observed GRB counts, M°^t Z2 , to be 
representative of the expected numbers, Af z xp _ z , we find 



\P*/Z1 — Z2 \ro, 



a/^* J>t£tM«0 (i + *) a 



M£\ J z *dz^(l + z)<* 



(4) 



7 In Kistler et al. 12008), where the focus was on intermediate z, a lower 
Li so threshold was used. 
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written in terms of the data for the groups of GRBs. The de- 
crease of (dV / dz) / (1 + z) at z > 1.5 increasingly amplifies 
the significance of distant observed GRBs. 

3. THE INFERRED HIGH-Z STAR FORMATION RATE 

We show our new determinations of the high-z SFR in 
Fig. Q] and how they compare with most of the conventional 
SFR results in this range 8 . While evolution with a = 3 (0) 
would decrease (increase) these values by a factor ~ 2, it 
would also disagree with z < 4 data. Other changes in the 
details of our analysis, while changing the results somewhat, 
do not change the important point that the high-z SFR must 
be relatively large to accommodate the several observed high- 
z GRBs. These details include the exclusion of a particular 
GRB or changes in redshift assignments or ranges. Taking 
into account the Poisson confidence interval for 4 observed 
events, and slightly enlarging it, we assign a statistical un- 
certainty of a factor 2 up or down. As the number of GRBs 
increases, this uncertainty will correspondingly decrease, and 
more attention will need to be placed on systematics 910 . Our 
luminosity cut happens to restrict most of the GRB redshifts 
used to have been obtained by absorption lines, reducing one 
possible source of error. 

We provide an update of the SFH fit of iHopkins & Beacoml 
(2006) based on our GRB results (solid line in Fig. |TJ. Here 
we adopt a continuous form of a broken power law, 



P*{z)=f>o 



(1 + z) a " 



1 + z 
B 



br/ 



i 



c 



(5) 



where a = 3.4 (b = —0.3, c = —3.5) is the logarith- 
mic slope of the first (middle, last) piece, and the normal- 
ization is po = 0.02 M® yr _1 Mpc -3 . Using r\ ~ -10 
smoothes the transitions, while r\ — » — oo would recover the 
kinkiness of the original form. The breaks at Z\ = 1 and 
z 2 = 4 correspond to B = (1 + z x ) 1 ^ a l h ~ 5000 and 

c = (i + zi)< b - a )/ c (i + z 2 y- b/c ~ 9. 

4. COMPARISON TO OTHER SFR RESULTS 

Our SFR results, particularly the z = 5 — 7 point, are at 
the high end of SFR results obtained with more conventional 
methods, which show significant differences among them- 
selves. Likely, these methods are measuring different pop- 
ulations of galaxies, with the GRBs probing low-luminosity 
galaxies that are otherwise not accounted for fully. As de- 
tailed, the statistics of the GRB data are adequate, and it is 
unlikely that our GRB results are overestimated, and if any- 
thing, the true SFR is even larger. We assume that a steep drop 
in the SFR is not being hidden by a steep rise in the fraction of 
stars that produce GRBs, beyond the evolution already taken 
into account; that would be even more interesting. 

The extent of extinction by dust on high-z SFR measure- 
ments is not yet strongly constrained. There are a number of 
indications that dust is ubiquitous over the range 4 < z < 6 

8 Thes e measur ements are based on the products of massive stars, and 
assume a Salpeter (1955) stellar initial mass functio n. A change in the low- 
mass IMF with z would necessitate a rescaling (e.g., Wilkins et al. 2008). 

9 While we use 8 z>4 GRBs from Butler et al. (2007), none were observed 
in the ~ 1 year since then, indicative of a possible systematic effect. This will 
not change our results more than the quoted uncertainties. 

10 Our normalizing to the whole range in z = 1—4 reduces the impact of 
a hypothesized deficit of GRBs in z = 1-2 (Bloom et al. 2003; Coward et al. 

l200alFioTe~et al.120071) 



(e.g., ICharv et al.ll2005t: lOuchi et al.ll2004l: lAndo et al.ll2006h . 
and the dust correction in this range can be up to a factor of 
~ 2 — 3. At higher redshift, there are no strong observational 
constr a ints. The dust corrections assumed in Bouwen s et al.l 
d2007i 120081) are generally small at high z; on the other 
hand, observations verify the ex istence of at least o ne heavily- 
obscured galaxy at z ~ 6.6 dCharv et al.l 12005), and UV- 
selected samples would be expected to be biased against such 
objects. We also note that the recent mid-IR detection of the 
progenitor of S N 2008S implie s that some short-lived stars do 
produce dust (P rieto et alj|2008t) . For all of the LBG-based 
SFR results in Fig. Q] we have included the dust corrections 
indicated by the respective authors. While the true dust cor- 
rections may be larger, it seems likely that an additional factor 
of several is implausible. 

The LBG surveys are most sensitive to the brightest galax- 
ies, and great efforts are made to define the UV galaxy lu- 
minosity function (LF) and how well it is sampled at the 
faint end. For examp le, the results shown in our Fig. [T] from 
iBouwens et al.l (120081) are integrated down 0.2 L r , for £* de- 
fined at Z = 3. For the three lowest-z points, they also in- 
tegrate down to 0.04 L*, yielding SFR results a factor of a 
few higher. Correcting these measurements to account for the 
contributions from even fainter galaxies is necessary but diffi- 
cult due to the poor observational constraints on the faint-end 
slope of the LF. Depending on the value of this slope, and the 
lower limit of the integration, the necessary correction may be 
as little a s a factor of tw o, or as much as an order of magni- 
tude (see lHopkinsll20061) . It also seems that Lya Emitters do 
not contribute significantly to the total SFR. 

5. DISCUSSION & CONCLUSIONS 

We developed an empirical method for estimating the high- 
z SFR using GRB counts, improving o n earlier estimates 
of the high-z SFR using GRB data (e.g.. iBerger et alj|2005 



Nataraian et all 120051; 



Jakobsson et 



[e^g Ben 
aU[2006; 



Le & Dermerl 



2007; Chary et al. 2007) in several ways, not least of which by 



using signifi cantly updated SF R ((Hopkins & Beacom 2006) 
and/or GRB JButler et alj2 007) data. Taking advantage of the 
improved knowledge of the SFR at intermediate z, we were 
able to move beyond the assumption of a simple one-to-one 
correspondence between the GRB rate and the SFR, account- 
ing for an increasing evolutionary trend. The higher statistics 
of the recent Swift GRB data allowed the use of luminosity 
cuts to fairly compare GRBs in the full z range, eliminat- 
ing the uncertainty of the unknown GRB luminosity function. 
By comparing the counts of GRBs at different z ranges, nor- 
malized to SFR data at intermediate z, we based our results 
squarely on data, eliminating the need for knowledge of the 
absolute fraction of stars that produce GRBs. 

Our GRB-based SFR value in z = 4 — 5 is comparable to 
more conventional results, which may be taken as a validation 
of our method, i.e., possi ble selection effe cts have been ade- 
quately controlled (see iKistler et al.ll2008l for more details), 
while our SFR value in z = 5 — 7 is larger, which has impor- 
tant implications. A significant population of low-luminosity 
star-forming galaxies picked out by GRBs but missed in other 
surveys could reconc ile these differences (searches utilizing 
cluster lensing, e.g., iRichard et alj 20081 may help, ifcos- 
mic variance can be understood). lYan & Windhorst! (2004) 
argue that the correction for the faint-end of the LF should 
be large to correct for dwarf galaxies missed in LBG surveys. 
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Sin ce GRBs are observed to favor subluminous host galax- 
ies dFvnbo et alj 120031: Ee Floc'h et al.1 120031; iFruchter et all 
2006), they may be probing such faint galaxies. 

The universe was f ully reionized by z ~ 6 dWvithe et al.l 
120051: iFan et al.l l2006h and it appears that AGN could not 
have been responsible (Hopkins et al. 2008). Can stars have 
reionized the universe? In Madau et al. (1999), the SFR den- 
sity required to produce a sufficient ionizing photon flux is 
parametrized by two factors: the photon escape fraction (/ csc ) 
and the dumpiness of the IGM (C). These enter as a ratio, 
C/fesc, which is not known precisely. For C// csc < 30, the 
required SFR at z = 6 is />» (z = 6) > 0.03 M Q yr" 1 Mpc -3 , 
just at the level of our GRB-based SFR. 

This ratio may be higher (due to either C or / esc ), how- 
ever, which would render our measured value of p*(z = 6) 
too low, reintroducing the concern that star formation may 
be insufficient to achieve reionization (e.g., Gnedin 2008, 
C. A. Faucher-Giguere et al. 2008, in preparation). Is there 
any way around this? One can plausibly increase our SFR de- 
termination by considering the increasing incompleteness of 
the GRB sample with z due to our i; so cut, since for fixed lu- 
minosity, hi gher- z bursts are relatively more difficult to detect 
dButler et al.l 120071) and the universe will at some z quickl y 
become opaque blueward of Lya (e.g.. lCiardi & Loebll2000l) . 
If there is a maximum metallicity for forming a GRB, then the 
evolutionary trend may saturate at high enough redshift; arbi- 
trarily assuming no evolution beyond z > 4 would give an ad- 
ditional factor of < 2. Finally, requiring a minimum amount of 
metals in a star for a successful GRB would suppress the GRB 
rate at increasingly high-z, necessitating a higher underlying 
SFR. A similar argument could be made concerning whether 
early, metal-poor stars possessed sufficient angular momen- 



tum. Overall, more data are needed concerning the varying 
ionization fraction with redshift to establish if star formation 
is quantitatively acceptable as an explanation for reionization, 
instead of simply by default if AGN have been eliminated. 

The possibility of continuing evolution in the GRB rate 
relative to the SFR invites some astrophysical speculations. 
While the overall black hole production rate in core-collapse 
supernovae is poorly understood empirically (Kochanek et al. 
2008), bright GRBs might be regarded as tracing the large 
angular momentum end of the black hole birth d istribution 
(e.g., for collapsars; iMacFadven & Wooslevlll999l) . Evolu- 
tion would mean that the high-z universe was more efficient 
at producing such black holes than usually considered. This 
may have implication s for the nucleosyn thetic yields from 
these explosions (e.g.. lNomoto et al]|2006l) . The effect of su- 
pernovae on the gas in high-z gala xies is an ongoing field of 
research (e.g., Whalen et al. 2008). Following from the con- 
siderations above, it may be that small galaxies had a dis- 
proportionately large rate of GRBs relative to normal super- 
novae. In this case, it would be interesting to examine the fate 
of such galaxies after including multiple injections of highly- 
asymmetric, relativistic ejecta, as well as implications for en- 
richment of the IGM (i.e., "GRB feedback"). 
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